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Abstract 

The radiative B — > frf, B — ► cuy decay modes are caused by the flavor-changing-neutral- 
current process, so they give us good insight towards probing the standard model in order to 
search for new physics. In this paper, we compute the branching ratio, direct CP asymmetry, 
and isospin breaking effects using the perturbative QCD approach within the standard model. 



1 1 Introduction 

> 



The standard model (SM) predicts large CP violation in B decays ^ |2] and they have been 



Cf) ' verified in B — > J/ipK s 011], B 7r7T El El; and B — > DK decays. The quest of high energy 

physics has always been to search for the most fundamental theory. So our immediate goal is to 
| search for deviation from the predictions of the SM. It is believed that the quantum effects in 

' B meson decay amplitudes may contain effects of new physics. 

The flavor-changing- neutral-current (FCNC) process which causes b — ► 37 and b — > dj decays 
Qh| may contain new physics (NP) effects through penguin amplitudes. As the SM effects represent 

Q_|' the background when we search for NP effects, we shall compute these effects. In doing so, we 

can understand the sensitivity of each NP search. 

The first experimental evidence of this FCNC transition process in B decay was observed 
about a decade ago, where the inclusive process b — > 57 and exclusive process B — > K*y 
^ ■ were detected, and their branching ratios were measured [S]. On the other hand, the expected 

branching ratio for b — > is suppressed by O(10 -2 ) with respect to that for b — > s-y, because 
of the Cabbibo-Kobayashi-Maskawa (CKM) quark-mixing matrix factor. The world average for 
b — > d penguin decays is given as follows 



Br(B° -> p°7) = (0.38 ± 0.18) x 10" 6 
Br(B° ury) = (0.541°;^) x 10~ 6 
Br{B + -> p+ 7 ) = (0.68±g;|?) x 10~ 6 . 



Theoretically, B —> and B — ► toy are widely studied both within and beyond the SM 
jl()l 111 j . The bound states are involved in the exclusive process, so the perturbation theory 
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Figure 1: The heavy b quark decays into the light d quark and a photon through the electro- 
magnetic operator, and the decay products dash away back-to-back with momenta 0(Mb/2). 
In order to form a p + meson with no hadron jets, the spectator quark must line up with d. This 
can be accomplished most efficiently by exchanging a hard gluon 



can not be used in a simple manner. It has been shown that, at least in the leading order, all 
nonperturbative effects can be included in the definition of the B meson and the vector meson 
wave functions, and the rest of the amplitude (the hard part of the amplitude) can be computed 
in the perturbation theory. This is called the perturbative QCD (pQCD) approach and it was 
proven several years ago j!2l I13j. In this paper, we compute the branching ratio, direct CP 
asymmetry, and isospin breaking effects for B — * fry, B — * ujj decays by using the pQCD within 
the SM. 

The remaining part of this paper is organized as follows. In Sec El we briefly review the 
pQCD approach, and in SecOJ we present some basic formulas such as the effective Hamiltonian 
and kinetic conventions. In SecEJ the hard amplitudes calculated in pQCD are given. Section 
[3] is devoted to numerical calculation and discussion. Finally, a brief summary is given in Sec® 

2 Perturbative QCD Approach 

In order to explain the pQCD approach, we want to suppose that a static B + meson decays into 
p + and 7 through the O-j^ operator as in Fig^ 

In the rest frame of the B + meson, the b quark is almost at rest and the spectator u 
quark moves around the b quark with O(A) = 0(Mb — va^) momentum, where Mb, and m& 
are B meson, and b quark mass, respectively. Then the b quark decays into d and 7, and these 
products dash away back-to-back with 0(Mb/2) momenta. When a quark is rapidly accelerated 
like this, infinitely many gluons are likely to be emitted by bremsstrahlung. There is a familiar 
phenomena in QED, when an electrically charged particle is accelerated, infinitely many photons 
are emitted. But the gluon emission by bremsstrahlung QCD must result in many hadrons in the 
final state. As the emitted gluon will hadronize, the fact that no hadron except for p(u) should 
be observed in B — * p(co)^, means that the bremsstrahlung gluon emission mentioned above 
can not occur. Thus the branching ratio for an exclusive decay B — * p(uj)j is proportional to 
the probability that no bremsstrahlung gluon is emitted. The amplitude for an exclusive decay 
contains the Sudakov factor and it is depicted in Fig|21 As seen in FigEl the Sudakov factor is 
large for small b and small Q, where b is the spacial distance between quark and antiquark into B 
meson, as shown in FigEl and Q is the b quark momentum inside the B meson. Large b implies 
that the quark and antiquark pair is separated in space, which in turn implies less color shielding. 
Similar absence of the color shielding occurs when the b quark carries the most of the momentum 
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Figure 2: The dependence of the Sudakov factor exp[— s(Q, b)] on Q and b where Q is the b 
quark momentum, and b is the interval between quarks which form hadrons. It is clear that the 
large b and Q region is suppressed. 




Figure 3: b is the transverse interval between the b and u quark in the B meson. 



of the B meson. That is, as seen in Fig|21 in order to form a p(co) meson with no hadron jets, 
the condition for color shielding is essential. The condition needed for the color shielding is the 
small separation in space between quark and antiquark within the meson, and it indicates that 
the energy scale of the decay process should be high. Actually, the invariant-mass square of the 
exchanged gluon depicted in Figgis about O(AMb), which can be considered to be in the short 
distance regime. Thus we can see that the decay process can be treated perturbatively. The 
decay amplitude for the exclusive mode like B — ► p(u)j decay can be factorized into the hard 
part with a hard gluon exchange, which can be treated perturbatively, and the soft part of all 
nonperturbative strong interactions is included in the meson wave functions. 
Then the total decay amplitude can be expressed as the convolution like 

/•l /-l/A 

/ dx\dx2 i d 2 b\d 2 b2 C(t) ® exp [— S(x\, x 2 , &i, b 2 , t)] 
Jo Jo 

®<&p,u>{ x 2,b 2 ) ®H(x 1 ,x 2 ,h,l>2,t)®&B(xi,bi), (!) 

where <&pw(x2,&2) anc ^ ^B^ii^i) are meson distribution amplitudes, exp [— S(x\, x 2 , b±, b2, t)} 
is the Sudakov factor, which results from summing up all the double logarithms of the soft 
divergences. H(xi,X2,b\,b2,t) is the hard kernel including finite piece of quantum correction, 
61, 62 are the conjugate variables to transverse momenta, and x\, X2 are the momentum fractions 
of spectator quarks. 

In the computation of the decay amplitudes with the pQCD approach, we adopt the model 
functions for the meson distribution amplitudes. The meson amplitudes are characterized by 
the strong interaction. The effective range of the strong interaction which can propagate, is 
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wide. Then the meson distribution amplitudes should be expressed as some averaged physical 
quantity. Thus the meson amplitude does not depend on the decay process etc. For the B 
meson wave function, we adopt a model For the p and u meson wave function, we use 

ones determined by the light-cone QCD sum rule ^Sj. The detailed expressions for the meson 
functions are in Appendix EI 



3 Basic formulas 

The flavor-changing b —* d'j transition induced by an effective Hamiltonian is given by 

+VM {c[ c \fj,)o{ c \») + C%\n)0f{n)} (2) 
-VtbV t * d { C i (^)O f (M) + C7 7 (Ax)0 77 (M) + C 8ff (^)0^(^)}]+(h.c.), 

i=3~6 

where Cj's are Wilson coefficients, and Oj's are local operators which are given by 
0[ 9) = (diqj)v-A(<ljbi)v-A, ( 2 ] = (diqi)v-A(q~jbj)v-A, 



Of' = (dibi) V - A ^2(qjqj)v-A, 0\ q ' = (dibj) v _ A ^(qj 



Qi V-A, 



<-[ 

5 } = (dibi)v-A^2{qjqj)v+A, Oq = {dib^v-A^iWijv+A, (3) 

7l = ^^^(1+75)6^, Sg = ^ m b diu^{l +76)2^6^, 

and we neglect the terms which are proportional to d quark mass in 77 and Os g . Here (qiqj)vTA 
means (^7^(1 z fy 5 )qj, and i, j are color indexes. With the effective Hamiltonian given above, 
the decay amplitude of B — > p(uj)^ can be expressed as 

A = (F\H eS \B) = ^J2Kb V idCi^)(F\O i ( l j,)\B), (4) 

*!<? 

where F denotes the final state p^ or ivy. In addition, the amplitude can be decomposed into 
scalar (M s ) and pseudo-scalar (M p ) components as 

A=(e* v - e;)M s + — ^— e^e^e*^ PPP V M P , (5) 

ry • r%y 

where Py, and P 7 are the momenta of p{to) meson, and photon, respectively, e* and e v are the 
relevant polarization vectors. The matrix element (F\Oi(p)\B) can be calculated in the pQCD 
approach. 

For convenience, we work in light-cone coordinate. Then the momentum is taken in the form 

P = [P ,P ,PT) = ^ -j= , -j= , [P ,P ) I , (6) 

and the scalar product of two arbitrary vectors A and B is A ■ B = A^B^ = (A + B~ + A~B + ) — 
A± ■ B±. In the B meson rest frame, the momentum of B meson is 

Mr 

P B = (P+,Pz,P B ±) = -|(1 > 1,0±), (7) 
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(a) > (b) 



C77 S u O7-, 



Figure 4: Contribution from operator 7l to P — * p(co)j decay. The photon is emitted through 
the operator, and hard gluon exchange is need to form a p{uj) meson. 

and by choosing the coordinate frame where the p or u meson moves in the "-" and photon in 
the "+" direction, the momenta of final state particles are 



Mb 
M B 



P 7 = (P+,P-,P 7± ) = -#(1,0, 0j_), (8) 



\ = (P+,Py,P v± ) = -^(0,l,0 ± ). (9) 



The momenta of the spectator quarks in B and p or to mesons are 



, _ -> Mb 
ki = (k{,k l ,k 1T ) = (-j=-xi,0,k 1T ), (10) 

+ - -> Mb 
k 2 = {kj,k 2 ,k 2T ) = (0, -j=x 2 ,k 2T ), (11) 

where x±, and x 2 are momentum fractions which are defined by x\ = kf /Pg, and x 2 = k^ /Py , 
respectively. 

4 Formulas of the hard amplitude 

In this section we give the amplitudes caused by each operator in Eq.®. 



4.1 Contribution of O 



n 



At first, we present the contribution of the electromagnetic operator 7l . The diagrams are 
shown in Fig0J In this case, the photon is emitted through the operator, and hard gluon 
exchange is needed to form a p(to) meson. Contributions of the 77 operator to the amplitudes 
M s and M p defined in Eq.© are as follows: 



M^ a) = -M^ a) 



x ry 



2P(°)& J d Xl dx 2 f d6 1 d6 2 6 1 6 2 a.(t?)exp[-5B(t?)-5v(t?)]5t(zi)C77 7 (t?)^B(a;i,&i) 
(K(x 2 ) + <P v (x 2 )\ H\ a) (A 7 b 2 , B 7 h,B 7 b 2 ), 

(t? = max(A 7> B 7 , l/6i,l/6 2 )), (12) 
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M 



S(b) 
7 7 



p(b) 



7-) 



-2F<% J d Xl dx 2 J d6id6 2 6i6 2 a s (t fe 7 )exp[-5B(4)-5y(t 6 7 )]^( a ; 2 )C 77 (4)0 B (x 1 ,6 1 

(l + x 2 )<t>l(x 2 ) + (1 - 2x 2 )r v [^(x 2 ) + 4> v v (x 2 )]\H\ h) {A 1 b ll C 1 b u C 7 b 2 ), 
t\ = max(A 7 ,C 7 ,l/bi,l/b 2 



H ( 7 a) (A 7 b 2 , B 7 b u B 7 b 2 ) = K (A 7 b 2 ) 6(h - b 2 )K (B 7 h) I (B 7 b 2 



-e^-b^KoiB^IoiB-rh) 



Hi b) (A 7 b u C 7 h,C 7 b 2 ) = H^{A 7 b u C r h,C 7 b 2 ), 
A 2 7 = x x x 2 Ml, B% = Xl M B 2 , C 7 2 = x 2 M B . 



(13) 



(14) 

(15) 
(16) 



Here Kq, Iq are modified Bessel functions which are extracted by the propagator integrations. 
We define the common factor as 



F (o) = %£c F Mj>, 
V2 vr 



(17) 



and the CKM matrix element as £ q = V* h V qs . The exponentials exp[— S's(t)] and exp[— Syit)] 
are the Sudakov factors [12], and the explicit expressions of the exponents Sb, Sy are shown in 
Appendix^! The quark structures for vector mesons are p + = \du), p° = \uu — dd)/v / 2, and 
oo = \uu + dd}/\/2, then the decay amplitudes for each decay modes caused by 77 operator are 
given as follows: 



M{BT 



p + i) 3 



M(B° -» A) J ' 



n 



1 

75 



i 



^ V2 

where j expresses the decay amplitude components S or P. 



(18) 

(19) 
(20) 



4.2 Contribution of 0$g 

The diagrams for the contribution of the chromomagnetic penguin operator Os g are shown in 
FigEI Contributions of each diagram are given in the following. In this case, a hard gluon 
is emitted through the Os g operator and glued to the spectator quark line, and a photon is 
emitted by the bremsstrahlung from the external quark lines. Each decay amplitude caused by 
Osq operator is expressed as follows: 



M 8 5(a) (Q 6 



-M 8 P(a) (Q b ) 



-F {0 ktQb J d Xl dx 2 J db x db 2 bxb 2 a s {tl)ex V [-SB{t a s ) - S v (t%)]S t (xx) 

x C% g (t%)4> B (x\, h) [x 2 r V (j)y(x 2 ) + xi(j)y{x 2 ) + x 2 r V (j) v (x 2 )] 
x H£\Azb 2 ,B s b x ,B%b 2 ), 

(ig = max(^ 8 , Bg,!/^, l/fta)), (21) 
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Figure 5: Diagrams for the contribution of the chromomagnetic operator 0% g . A hard gluon 
is emitted through the 8fl operator and glued to the spectator quark line. Thus a photon is 
emitted by the bremsstrahlung of the external quark lines. 

M* (b \Q d ) = -M^ b) (Q d ) 

= -F {0) ZtQ d J d Xl dx 2 J db^hhasit^eM-SBii) - S v (t b 8 )]S t (x 2 ) 

x C& g (tl)(f>B(xi,bi) [-3x 2 ry</>y(x 2 ) + (2x 2 - x 1 )(py(x2) - 3x 2 r V (py(x 2 )] 
x J fff ) (^ 8 & 1 ,C 8 6 1 ,C 8 6 2 ), 

(t| = max(^ 8 ,C 8 , 1/61, 1/62)) , (22) 

M! (c) (Q q ) = -M 8 P(C) (Q 9 ) 

= -F^ \tQ q J dxidx 2 J dhdhhha^expi-SBitl) - Sv(t c 8 )}S t ( Xl ) 

x Cs g {tl)(f>B{xi, h) [x 2 r V (py(x 2 ) - xi0y(x 2 ) + x 2 r v cpy(x 2 )] 

x H { s c) y\Ai\b 2 ,D s b u D s b^ 

^§ = max(^4^|,£) 8 , 1/61, 1/63)) , (23) 

Ms (d) (Q q ) = -F<-%Q q j ' d Xl dx 2 J dhdhb^a^eM-SBitt) - S v (t d 8 )]S t (x 2 ) 
x C S g(ti)<j> B {xi, 61) [(x 2 - xi + 2)<#(x 2 ) + 6x 2 r y ^(x 2 )] 
x Hf(j\Ai\b u E s b u E s b 2 \ : (24) 



M 8 (d) (Q q ) = FW&Qq J d Xl dx 2 J d6id6 2 6i6 2 a s (tg)exp[-5 fl (tg) - S v (t d 8 )]S t (x 2 ) 
x C Sg (4)^ B (xi,bi) [{x 2 -xi + 2)<$(a? 2 ) + 6x 2 r V (/)^(x 2 )] 



x H t 



\A^\b 1 ,E s bi,E a b 2 

(tj = max(^4^,£ 8 , l/b u I/62)) , 



(25) 



H ( 8 a) (A 8 b 2 ,B 8 b u B 8 b 2 ) = K (A 8 b 2 ) 9{bx - b 2 )K (B 8 bx)I (B 8 b 2 ) + (6 X 63) , (26) 



H^{A 8 b u C 8 b u C 8 b 2 



K Q (A 8 h) 9(h - 6 2 )fW(C 8 6 1 )J (C 8 6 2 ) + (6 X <-> 6 2 ) , (27) 



^ c) (J|A«|6 2 , J D 8 6i,D 8 &2) = ^K 2 ) #o( \l\A'i\b 2 ) eih-b^KoiDsbMD^ + ih^h 



+9(-A'i) i|^ 1) (^/i^|6 2 ) [»(6i - b 2 )K Q (D 8 bi)I (D 8 b 2 ) + (61 <-» 6 2 )] , (28) 



F^ dJ (^/|^|6i,^6i,^62) = ^4 2 ) *|ifo(^|^|6i) - 6 2 )^ 1) (^ 8 6 1 )J ( J E 8 6 2 ) + (6 X <-> 6 2 ) 



0(-Af) l~Y Hg-\j\A*\h) dih-b^H^iE^ME^ + ih^b,) 



C 1 )/ 



(29) 



= xix 2 M|, S| = M|(l + xi), Cf = M|(l-x 2 ), 
" 2 = (x 1 -x 2 )M|, 



.4 



L»f = Xl M B 2 , El = x 2 M B 2 . 



(30) 



Here we define Q q as the electric charge for the quark q: = Q d = — 1/3 and Q u = 2/3. Then 
the decay amplitudes for each decay channels can be written as follows: 



M{B + -» p+ 7 ) J 



M^ a) (g fc ) + M^>(Q d ) + M^(Q U ) + M|W(Q tt ), 



(31) 



M(B° -» p° 7 ) j 



1 

71 



1 



^2 



'8s 



(32) 
(33) 



4.3 Loop contributions 

In this section, we consider the contributions of diagrams with the effective operators Oj's 
inserted in the loop diagram. 0\ does not contribute because of the color mismatch. Penguin 
operators 03^6 insertion is neglected, because they are small compared with 2 insertion in the 
loop diagram. Therefore, we only consider the tree 2 operator insertion. These diagrams can 
be separated into two types. One type is that of a photon emitted from the external quark line 
(FigEJ), and the other is that of a photon emitted from the loop quark line (FigJTJ). 



8 




Figure 6: Diagrams in which the operator 2 is inserted in the loop, and a photon is emitted 
from the external quark line. 0\ does not contribute and it can be shown that 03^6 can be 
neglected. 



4.3.1 Contributions of external-quark-line emission 

For the calculation of the diagrams in FigEl one can at first calculate the effective vertex b — > dg 
by performing the loop integration. For the topological structure with 2 inserted in the loop 
diagram of FigEl the effective vertex obtained with MS scheme is 



V 



~-G(mlk 2 ,v) 



dx4x(l — a;) log 



mf — x(l — x)k 2 



te 



(34) 
(35) 



where i = u, c is the flavor of the loop quark, k is the momentum of the virtual gluon, and v is 
the Lorentz index of the gluon field. We can see that the vertex function has gauge invariant 
form. With the effective vertex given in Eq . (|34j) . the contributions of diagrams in FigHOcan be 
obtained as follows: 



M^ a) (Q b ) = M? ia \Q h 



\F^iiQ b I d Xl dx 2 I dhdbihhasitDeM-Mtt) ~ Sv(t%))S t ( Xl ) 



x C 2 (t a 8 )Mxi,bi)x 1 x 2 r v [^ v (x 2 )-^(x 2 )}Hi a) (A 8 b 2 ,B 8 b 1 ,B 8 b 2 ) 

2' 

G(m 2 , -xix 2 m 2 B ,tl) - - 



(36) 
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Mf\Q d ) 



l;F^CiQ d [ d Xl dx 2 [ db 1 db 2 bib 2 a s (tl)eM-SB(4) - S v (t b 8 )}St(x 2 ) 



xCtfDfofaM) [3x lX2 <l%(x 2 ) + x 2 2 r v {cj> v v (x 2 ) + cj>^(x 2 )}] flf^&i, C 8 h,C s b 2 ) 

(37) 



G{mi,-x lX2 m 2 B ,t^ 



M s ^ c \Q q ) 



= \F^iiQ q f d Xl dx 2 f dhdhhha^eM-SBitt) - S v (t c 8 )]S t ( Xl ) 

x C 2 (f % )<f> B {xxM) [x2T V {<Pv(x2) +<f>v(*2)} -x 1( j%(x 2 )] (J\Ag\h, D 8 b u D 8 b 2 ) 



G (m i ,( X2 - xi)m B ,tl) 



(38) 



lF {0 kiQ q [ d Xl d X2 f dbidhhhasit&expi-SBiti) - Sv(4)]S t (x 2 ) 



2' 

x C 2 (ti)0B(zi, h) [3(x 2 - xi)^(x 2 ) + X2 r v {3(l + x 2 )<t> v v {x 2 ) - (1 - x 2 )^{x 2 )}] 

G(ml(x 2 - Xl )rn 2 B ,tt)-~ , (39) 



x H^UlA'Ub^Esb^Esb,) 



Mu (d) (Q q ) = ~F<®SiQ q j ' dx x dx 2 j dhdhhha^tDexpi-SBiti) - S v (t d 8 ))S t (x 2 ) 



x C 2 (4)</> B (xi, h) [3(x 2 - xx)<t^{x 2 ) - x 2 r v {(l - x 2 )4> v v {x 2 ) - 3(1 + x 2 )ffi(x 2 )}] 

2' 
3 



x Hi d) U/\A' 8 2 \b 1 ,E s b u E 8 b 2 



G (m 2 , {x 2 - Xl )m 2 B ,tfj 



(40) 



Then the decay amplitudes in this case can be expressed as follows: 



M(B + -> p+7) J 



ii 



M(B° -> p° 7 ) j 



1/ 



1 

72 L 



M^ (a) (Q fc ) + M/f } (Q d ) + M^ (c) (Q d ) + M^(Q d ) 



1 



(41) 

(42) 
(43) 



M(B° - ^ = -±= [Jlf (Q 6 ) + Mjf (Q d ) + M x f ^ (Q d ) + M$*> (Q d ) 
4.3.2 Contributions of internal-loop-quark-line emission 

The diagrams in which a photon emitted from the internal loop quark line are shown in Fig0 
The sum of the effective vertex b — > djg* in FigsJTfa) and[7Jb) has been derived in |17l I18j . The 
result can be expressed as 

I = - 75 )T a bJ w ^ £ ;, (44) 

with the tensor structure given by 



■ A A [(q ■ k)e fiupa (q - k) a + e upar q a k T k fM - e ppar q a k T q u ] 

+^4-5 [^iiparQ k k v — k Gfiupo-q ] , 



(45) 
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(a) 



(b) 





Figure 7: Diagrams in which the operator 2 is inserted in the loop, and a photon is emitted 
from the internal loop quark line. 



and 



A, 
A 5 



xy 



x(l — x)k 2 + 2xyq ■ k — mf + ie 
x(l — x) 



x(l — x)k 2 + 2xyq ■ k — mf + ie 



(46) 
(47) 



where q is the momentum of the photon q = P B — Pv, arid k is the momentum of the gluon 
k = k 2 — k\. The result of the amplitudes M s and M p contributed by each diagram in Figd 
can be expressed as 



Mi 



.l_p(o)£. J dx J dy J dxidx 2 j dbibia s (t 2i ) exp[-S B (t 2i ] 



x C 2 (t 2i )4> B {x l M)H2i{b l AM\fW\)- 



xyx 2 M B — mf 
x(l - x)x 2 (3x 1 (py(x2) + x 2 r v {4>v(x 2 ) + 0y(x 2 )}) 
-xyx 2 ((1 + 2xi)<$(x 2 ) - r v {(l - 2x 2 )<j) v v (x 2 ) + <fy(x 2 )}) 



(48) 



M 2 P t 



1 rl—x r r 

dx dy dx\dx2 \ d&i&ia s (i 2i ) exp[-Ss(t 2i )] 
a Jo J J 

1 



x C2(t2i)Mxub 1 )H 2i (b 1 A,b lv W\)- 



xyx2Mg — mf 
x(l - x)x 2 (3xi^(x 2 ) + X 2 r v {(j)\ / (X2) + 0y(x 2 )}) 
-xyx 2 ((1 + 2x x )^{x 2 ) - rv{(l - 2x 2 )4> a v {x 2 ) + <fy(x 2 )}) 
(t 2i = max(Av^|,V&l)) , 



(49) 



A 2 = xix 2 M B , 



B 2 = xix 2 M B 



1 — x 



x 2 M 2 B + 



/??7 



x(l 



(50) 



H 2i {bxAMy/W\) = K (hA) - K (b lv W\) 
= KoihA)-^ 



(B 2 > 0), 



H (h y/\B2\) (B 2 <0). 



(51) 
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(c) (d) 




Figure 8: Annihilation diagrams in which the operators 0±, 2 are inserted. The box denotes 
operator insertion. 



and the decay amplitudes for each decay modes are as follows: 

M(B+ -> p +i%. = M 3 2i , M(B° - AH, = -^f , M(B° -> ^ = (52) 
4.4 Annihilation diagram contributions 

Next we consider the annihilation- type diagrams. They provide the main contribution for the 
isospin breaking effects in Br(B + —* p + ^f) and 2Br(B° —* p 0r y). 



4.4.1 Tree annihilation 

We consider the tree annihilation caused by 0%, 2 operators shown in Fig|HJ 

In the charged mode, this contribution is color allowed; on the other hand, it is color sup- 
pressed in the neutral modes. We define the combinations of the Wilson coefficients as 

oi(t) = Ci(t) + C 2 (i)/3, a 2 (t) = C 2 (t) + Ci(t)/3, (53) 

and each decay amplitudes can be given as follows: 



M S £\Q b ) = M^\Q b 

3y/6Q b f V TT 



-F {0) tu ®M* V * TV I dXl j dh ^MSB(t a A )]S t { Xl ) 
x o*(ti)0B(zi,6i)X'o(6iA.), (ti = max(A ,l/6i)), (54) 



Mf } (Q 9 



4M| 



dX2 



d& 2 6 2 exp[-S , y(t^)]5 4 (x 2 ) 



ajfc(t 



|ff«(6 2J B a ) [x 2 <^(x 2 ) + (2 - x 2 )^(x 2 ) 



(55) 
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M P A ^(Q q ) = F^ku 3 ^^ r v J dx 2 j db 2 b 2 eM-Sv(t b A )]S t (x 2 ) 

x a k (t b A )i^H^\b 2 B a ) [(2 - x 2 )^(x 2 ) + x 2 <fy{x 2 )] , 
t\ = max(B a , l/b 2 



(56) 



AMI 



dxi J dbibi exp[-S B (t A )]S t (xi) 
x akitftifoixuh^oihCi), (t° A = max(C a , , 



(57) 



M S A f{Q qll ) 



= F(°ku 3 ^'/ B7T rv J dx 2 J db 2 b 2 eMSv(t d A )}S t (x 2 ) 

x a k (t d A )i^H^(b 2 D a ) [-(1 - x 2 )<ftr{x 2 ) + (1 + x 2 )(f> v v {x 2 )\ , (58) 



< fc W (Q 9 ") = F {0) Cu ^/ B7r r v J dx 2 J db 2 b 2 exp[-S v (t d A )]St(x 2 ) 
x a k (t d A )i^H^\b 2 D a ) [(1 + x 2 )<ftr(x 2 ) - (1 - x 2 )^(x 2 )] , 



(4 = max( J D a ,l/6 2 )) , 



(59) 



A 2 a = (l + Xl )Ml 



B 2 a = {\- x 2 )M 2 B , 



C 2 a = Xl Ml 



D 2 a = x 2 M%. (60) 



Here we use the index k in order to express the Wilson coefficient combination in Eg . (|53j) . Then 
each decay amplitude can be expressed as follows: 



M(B~ 



P + l)\ = Mi a \Q b )+M A f(Q d )+M£\Q u )+M A ^(Q u ), 



(61) 



^ _ y M^\Q b ) + M^f (Q u ) + M^\Q d ) + Mlf (Q u ) 
M(B° - u 7 y A = ± [M j }^ (Q b ) + M*® (Q u ) + Mlf (Q d ) + (Q u ; 



(62) 
(63) 



4.4.2 QCD penguin annihilation 

Next we consider the QCD penguin annihilation contributions. There are two types of the 
annihilation diagrams in which the operators Oj's are inserted. One type is shown in FigEOand 
the other is in Fig llUI 

First we consider the FigOU contributions. Here we also define the combinations of the Wilson 
coefficients as 



a 3 (t) = C 3 (t) + C 4 (t)/3, 
a 5 (t) = C 5 (t) + C 6 (t)/3, 



at(t) = C 4 (t) + C 3 {t)/3, 
a 6 (t) = C 6 (t) + C 5 (t)/3. 



(64) 



For (V — A)(V — A) operators, the results are as follows: 
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(a) 




Oi 







Figure 9: Annihilation diagrams in which the QCD penguin operators are inserted. The box 
denotes operator insertion. 



M^f (Q b ) = M%? 

= F(0) ^ 3 ^2 /y V / / dhheMSB^St^) 
x a k {t a Al )(t)B{xiM)Ko{biA a ), (t A1 = max(i a , I/61)) 



(65) 



(o s ) 



d&afcexphM^i)]^) 



x a k (t b Al )i^-H^(b 2 B a ) [x 2 <jfy(x 2 ) + (2 - x 2 )^(x 2 )] , 



(66) 



M 



p(py 



(Q q ) = -F<>% 



SVGQqfBK 



4M| 



-TV 



J da* J 



db 2 b 2 exp[-S v (t b A1 )]S t (x 2 ) 



x a^)^ 1 ^^) [(2 - X2)Mx 2 )+X2Mx 2 )] , 



(t b A1 = max(£ a , l/6 2 )) , 



(67) 



F^ ^^^ rv J d Xl J dhheMSB^Stix,) 



4M| 



x a fc (i^ 1 )(Ma;i,6i)li'o(&iCi), (t^i = max (Ca, l/^i)) , 



(68) 



<{f (Q q >) 



3\/6Q q/ f B 7r 



4M| 



d6 2 & 2 exp[-Sy(4i)]StM 



x a k {t d A1 )i^H^\b 2 D a ) [-(1 - x 2 )0 y (x 2 ) + (1 + x 2 )4> v v {x 2 )] 



(69) 
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m2?~(Q<) = -F {0) ^^/ B7r r v J dx 2 J db 2 b 2 eMSv{t d A1 )]S t ( X2 ) 



x a k (t d A1 )i^H^\b 2 D a ) [(1 + x 2 )<t>\r{x 2 ) - (1 - x 2 )4> v v {x 2 )] 

*li = max(D , l/6 2 )), (70) 



A^ = (l + xi)Mi, Bi = (l-x 2 )M^, Ci = x x M^ B l a = x 2 M l B . (71) 

Here we use the index in order to express the Wilson coefficient combination in Eq.(|64|). and 
upper index "-" means the (V — A)(V — A) vertex structure. The total amplitudes in FigUJ] 
contribute only to the neutral decay modes and they are color suppressed contributions, thus 
they are given as follows: 



1 



M(B° -> urj) 



Al 



V2 
1 

72 



13 



ri(6)- 
13 



\Q d j) + (i<f (Q B ) - (Q d )) , (72) 



2Mif (Q 6 ) + {Miff (Q„) + Miff (Q d )} 

+2M*jf (Q d ) + {i\4f (Q u ) + i<f (Q d )} ] . 



(73) 



Amplitudes with (V — A)(V + A) operators can be related to those with (V — A)(V — A) 
amplitudes as 





+ (Q b ) 


= M^ a) ~ 


(Qb), 




(Qb) = 




(74) 




\Q q ) 


= 


(Q q ), 


mT 


(Q q ) = 


-M^'(Q q ), 


(75) 


M S A f 


(Q g ') 


= M S A ^ 


(Q q >), 




\Q q >)-- 




(76) 




iQg") 


= M S A f Y 


(Qq")i 




+ (Q q ») 


= -M^ (Q q „), 


(77) 



where "+" expresses the (V — A)(V + A) vertex structure. The decay amplitudes caused by the 
(V — A)(V + A) vertex exist only in the neutral modes and can be expressed as follows: 

M(B° -> A)£ = j= [{Mff (Q u ) - M^(Q d )} + {M^\ Qu ) - M A [f (Q d )}] ,(78) 



M(B° — ► wy)^ = ^[2M A [f (Q b ) + {M^ + (Q u ) + M^f (Q d )} 



-2M^ + (Q d ) + {i<f (Q u ) + M A [f(Q d )} 



(79) 



Next we consider the type two diagrams shown in Fig llOl For (V — A)(V — A) operators 



inserted in these diagrams, the results M A2 are the same as M A ^ 



M(B + -> p+7) 



.42 



M^l (Q b ) + M^l' (Q d ) + M*J£- (Q u ) + MM" (Q u 



L A2 4 



A2 A 



L A2 4 



(80) 



M (B» - p%Y A2 = ~ [M^l (Q b ) + M'Z[ (Q d ) + (Q d ) + (Q d )\ , (81) 

1 



ri(c)" 



M(B» -+ u 7 y A2 



V2l 



i<f (Q b ) + MZi (Q d ) + M A % (Q d ) + Mffi (Q d ) . (82) 



2 4 



3(c)- 

A 



2 4 
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Figure 10: The other type of annihilation diagrams with operator insertion. 

On the other hand, for (V — A)(V + A) operators, the results are 

M S A ^ + (Q b ) = M^f (Q h ) = 0, 



M 



S(6)H 



(Q d ) = -M S A ^ + (Q d ) 



F (0) e 



p 2M| 



dx 2 / db 2 b 2 exp[- S v (t A1 )]S t (x 2 ) 



x a k {t h A1 )4>l{x 2 )i^H^\b 2 B a ), 



M S A £ + (Q q ) = M^ + (Q q ) = 0, 



M s A t(Q q 



-M 



s(dr 



A2 k 



(Qa 



= F^ ^^I^ I dx 2 I db 2 b 2 eM-Sv(t d A1 )]St(x 2 ) 



2M| 



x a k (t d A1 )^(x 2 )i^H^\b 2 D a ), 
then the amplitudes of this type with (V — A)(V + A) vertex become as follows: 
M(5+ - „+ 7 )£ = M A [f(Q b ) + M^ + {Q d ) + Mlg + (Q U ) + M™ + {Q U ), 



V2 



i<f (Q b ) + M'Zl ' (Q d ) + M'Jg ■ (Q d ) + ' (Q d 



f 3(b) 



j(cy 



M(B° - ^)\\ = j= [i<f (Q b ) + (Q d ) + Mjg + (Q d ) + M^f (Qd) 

In these type two cases, they are all color allowed decay modes. 
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4.5 The final decay amplitudes M s and M p 



Finally, we summarize the amplitudes M v s, j = S,P for each decay mode: 
M j (B + ^p+j) 



M(B + - 




- 1 y 7l + M{B + 




; + M(B+ 


- P + l) 3 u 


+M(B + - 




- 1 y 2i + M(B+- 




+ M(£+ - 




+M(B + - 


^p" 


' 









(90) 

M^(S°^A) = M(B° — ► /J°7)7 7 + M(B° — ► p°7)i g + M(B° — ► p°7)^ 
+M(5° - A)4 + M(B° - p° 7 )^ + M(5° - p°7)^ 
+M(5° - p° 7 )£ + M(£° -> p° 7 )^ + M(£° - p°7)^ 2 , (91) 

JVP(B — > ury) = M(5° -► u; 7 ){ 7 + M(5° w 7 ) J 8g + M(5° u; 7 )^ 
+M(5° -► io^y 2i + M(£° -► w 7 )^ + M(B° -► u; 7 )^ 
+M(5° -» w 7 )£ + M(5° -> W7 )i" 2 + M(S° -» W7)^ . (92) 



5 Numerical analysis and discussions 

In our numerical calculations, the choice of the input parameters is summarized in TabQ where 
A, A, p and r\ are CKM parameters in Wolfenstein parametrization ^j, and p = p(l — \\ 2 ), 
fj = r](l — |A 2 ). Their values can be found in PDG |20| . The numerical results for each decay 
amplitudes M? in the B° -» p°~f (TabEJ), B° -» loj (TabEJ), and B + -> p + -/ (TabgJ) in unit of 
lCT 6 GeV~ 2 are as follows. 

When we estimate the physical quantities like branching ratio, direct CP asymmetry, and 
isospin breaking effect, we take into account the following theoretical errors. The detailed 
discussions for the errors are in |21| . First, we change the input parameters; the decay constants 
and ujb in. the B meson wave function, and we regard the 15% error in each cases at the amplitude 
level. This generates the theoretical error for the physical quantities about 40% in the branching 
ratio, 5% in the direct CP asymmetry, and 30% in the isospin breaking. 

Second, we estimate that the higher order effects in perturbation expansion to be about 
15% error in the amplitude. This leads to about 30% in the branching ratio and in the isospin 
breaking. Here, the cancellation of the higher order effects can occur by taking the ratio of the 
decay width in the direct CP asymmetry. Then we can neglect these uncertainties for the CP 
asymmetry. 

Third, the error due to the CKM parameter uncertainties p and fj generates about 30% 
error in the branching ratios and direct CP asymmetry, and 100% error in the isospin breaking 
effects. We can see that the uncertainty which comes from the CKM parameters are large 
compared to B — * K*j decay modes (2^. The reason for it is that the all CKM matrix elements 
which concern B — > p{oj)^ decay (V^Vtd, V* b Vcd, V* b V u d) are comparable, and the three angles 
of the KM unitarity triangle are sizable: the situation is different from in B — > K*j decay. The 
conditions mentioned above make the uncertainty from the CKM parameters large. 

In the end, we also take into account the uncertainties from u quark loop contributions like 
FigsEJandJTJ We guess that the nonperturbative effects in the u quark loop might lead to large 
hadronic uncertainties. So, we introduce the 100% theoretical error at the amplitude level. This 
theoretical uncertainty leads to small uncertainties (about 2%) for the branching ratio, 80% 
errors for the CP asymmetry, and about 3% errors for the isospin breaking effects. 
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A^ 4) 

MS 



CKM parameters and QCD constant 
A A p fj 

0.2196 0.819 0.20 ±0.09 0.33 + 0.05 250 MeV 

Meson decay constants 

fp fp fu) foj 

220 MeV 160 MeV 195 MeV 160 MeV 



fB 

190 MeV 



M w 
80.41 GeV 



M B 
5.28 GeV 



Masses 
M p 
0.77 GeV 



0.78 GeV 



1.2 MeV 



B meson life time 

T B T B ± 

1.542 ps 1.674 ps 



Table 1: Summary of input parameters 



















-172.12 


-0.44-1. 19i 


-7.76 - 3.45 i 


172.12 


0.48+1.181 


7.55+ 3.47 i 


M» c /F^i c 












0.39+l.Oli 


-1.21+8.84i 




-0.08-0.97i 


0.42-6. 281 




M?JF^ U 












1.35+2.06i 


-l.ll-27.76i 


1.14-O.Oli 


-1.13-1.98i 


-0.73+28.15i 


-2.47+0.17i 


Table 2: The numerical results for B° - 


->■ p°7 decay at p = 0.20, rj = 


0.33, lo b = 0.40GeV. 


















161.59 


0.44+1. 21i 


7.73 + 3.45 i 


-161.59 


-0.47-1. 181 


-7.71 - 3.42 i 














-0.46-0.99i 


1.21-8.581 




0.12+0.99i 


-0.36+6.35i 




Ml/F^Cu 












-1.03-2.10i 


1.64+ 27.00i 


1.04-O.Oli 


1.06+2.10i 


0.05-27.45i 


-2.25+0.12i 


Table 3: The numerical results for B° 


— > 0^7 decay at p = 0.20, r/ = 


0.33, to b = 0.40GeV. 
















Zi^MZjFW^ 


243.72 


4.76-3.121 


-4.61 - 2.73 i 


-243.72 


-4.56+3.15i 


4.08 + 2.42 i 




W c 






W« 




-0.74+2.701 


1.70-15.361 




1.51-3.021 


-0.48+ 11.231 
















-2.39+ 6.05i 


1.70+39.28i 


37.28-7.90i 


2.65-5.191 


0.99-39.80i 


-55.47-0.37i 



Table 4: The numerical results for B + -> p+7 decay at p = 0.20, fj = 0.33, w B = 0.40GeV. 
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Then the total theoretical error for each physical quantity becomes as about 60% in the 
branching ratio, 85% in the CP asymmetry, and 100% in the isospin breaking effects. 

With the amplitudes M and M p defined in Eq.Q, the total decay rate of B — > p(u>) 7 is 
given by 

1 " 8nM B ' ^ 
and the relevant decay branching ratio is defined to be 

Br = , (94) 

where tb is the mean lifetime of the B meson. The branching ratios for neutral and charged 
modes are defined as 

Br(B ± -» p ± 7) = i [Br(B + -» p+ 7 ) + Br{B~ -» p~ 7 )] , (95) 
Sr(S°^p 7 ) = ^ [Br(B° -» p° 7 ) + 5r(5° -» p° 7 )] , (96) 

and its' predicted values become as 

Br(B° -»■ p° 7 ) = (1.2 ±0.7) x 10" 6 , (97) 
Br(£° -^ury) = (1.1 ± 0.6) x 10~ 6 , (98) 
Br(B ± p ± 7 ) = (2.5 ± 1.5) x 10~ 6 . (99) 

The direct CP asymmetry is defined by 

± + y) _ r(B- -» p- 7 ) - r(g+ -» p+ 7 ) 
^ "* ' 7) - r ( 5-^p- 7) + r(5 + ^p +7 ) (100) 



for charged B meson decays, and 
A cp (B° - p» 7 ) 



o _ fVv ^ = r(B° -> P°H 7 ) - T(E° -» p_Vh) 

r(B° -► pVb) + r^B -► p°(w) 7 ) 



Iioi) 



for neutral -B meson decays. The numerical results for these CP asymmetries in B — > p 7 and 
w 7 are as follows: 

A cp (B° -» p° 7 ) = (17.6 ± 15.0)% (102) 



4:p(-B -» w 7 ) = (17.9 ± 15.2)% (103) 
Ay^-^A) = (17.7 ± 15.0)%. (104) 

Next we discuss the isospin breaking effect in B —* p 7 decay. The isospin relation requires 
that the branching ratio of B + — ► p + 7 is two times of U° — > p° 7 . However, the contribution of 
the annihilation diagrams can violate this isospin relation. We can define the isospin breaking 
parameter as 

Ao + (5-P7) = #77^4 "I. (105) 

= 2 r r^oZy 7 ) ) w 

A(p 7 ) = A ° + : A °- . (107) 
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Numerical results 





J3IdllCJ.lJ.Ilg IdllO 




B° -» flS 


R0 ^ 

i J r to 1 y 


p j 


fi o _i_ n 71 y If) -6 

I J. . _L_ VJ. 1 J A -L L/ 


t-\ i _|_ n a\ v i n— 6 

I ± . ± _l_ U.UI A -LVJ 


(n K i i K\ v i rv— 6 

li.tj _1_ J..U! A ±VJ 




T^irpff asvininptTV 

L/ 11 VjV^ \j V_y X CAjO y llllllv tl y 










(17.6 ± 15.0)% 


(17.9 ± 15.2)% 


(17.7 ±15.0)% 




Isospin breaking effects 






A(p 7 ) = -(5.4 ±5.4)% 





Table 5: The conclusion related to the branching ratio, CP asymmetry, and isospin breaking 
effects. 



If isospin relation is maintained, AQ97) defined above should be zero. Our numerical result for 
isospin effects is 

A(p7) = -(5.4 ±5.4)%. (108) 

6 Conclusion 

In this paper, we calculated the branching ratio, direct CP asymmetry and isospin breaking 
effects within the standard model using the pQCD approach. Our predictions for the physical 
quantities are summarized in TabJSJ 

The B —* p magnetic form factor T p is defined as 

(p(P2,e K *) I U? da p,b I B(P X )) = -iT»($)e m p p e a p pPq p (109) 

where P = P\ + P2, q = P\ — P2, and the value computed by the pQCD approach is = 
0.26 ± 0.07. The value from the light-cone-QCD sum rule (LCSR) is T( = 0.29 ± 0.04 22 , then 
our value of the form factor is in good agreement with LCSR. The branching ratios only from 
77 become Br(B ± -► ^7) = (2.4 ± 1.2) x 10~ 6 , Br(B° -► p°-f) = (1.1 ± 0.5) x 10~ 6 , and 
Br(B° —* loj) = (1.0 ± 0.5) x 10~ 6 ; then by comparing them to Eqs.l(57 jl -(|55 jl . we can see that 
O-j-f contributions are dominant. 

The subtle excess of the branching ratio of B° — > /?° 7 compared to that of B° — ► loj is caused 
by the following two reasons: (1) the difference in the meson mass and decay constants between 
p and lo; (2) the annihilation contributions from 0\ to Oq. We examined these possibilities, and 
concluded that the subtle excess of the branching ratio mainly comes from (1), and the effects 
from (2) are very small. 

The isospin breaking effect A(pj) is caused by the contributions Os g (FigEJ, charm and 
up quark loop contributions (FigEJ), and 0\ ~ Oq annihilation contributions (Figs l8llU|) . Our 
prediction for this quantity is given in Eq. H108j) . The main contributions to the isospin breaking 
effect come from Oi ~ Oq annihilation diagrams. In general, we can expect that the annihilation 
contributions are suppressed by the factor 0(m q /mb), where m q = m u , m^. In our computation, 
weak annihilations caused by O3 ~ Oq are about 5% and tree annihilations caused by 0\ , O2 are 
1% in the neutral modes, (see Tabs |2] and EJ); on the other hand, in the charged mode, weak an- 
nihilations are about 2% and tree annihilations are about 20%; this contribution is large because 
it is a color allowed process (see Tab0J) in the amplitudes. If we neglect 0\ ~ Oq annihilation 
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contributions, the isospin breaking effects have the opposite sign: A(/ry) = +(3.4 ± 3.4)%. Thus 
the annihilation contributions are crucial to the isospin breaking effects. 

When we compare our results with the world averages of experimental data for the b — > dj 
decay modes §||, our results for the branching ratios are somewhat large. For now, we shall not 
worry about it for the following reason: Note that our conclusion given in TabJHl 

Br(B° p°j) ps Br(B° ury) ~ ^Br(B + p+j), (110) 

follows from the isospin symmetry and the fact that contribution from the Oj^ operator dom- 
inates over all other contributions. A similar conclusion has been derived from the B — > K*^ 
decay mode |23| . and experimental results for B — > K*^ agree with our conclusions. While the 
error is large, the relationships indicated by Eq. (jllO|) are not obviously seen in the recent ex- 
perimental data. We thus feel it is too early to discuss the validity of Ea. ()llUj) , We expect that 
the data may change by about a factor two if Eq. ljllOj) is approximately valid. 
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A Wave function 

In the calculation of the decay amplitude, the wave functions of meson states can be defined 
via nonlocal matrix elements of quark operators sandwiched between meson states and vacuum. 
Next let us introduce the wave functions needed in this work. 

The two leading-twist B meson wave functions can be defined through the following nonlocal 
matrix element |24| : 



d 4 z 



iki z 



o (2vr) 4 
i 



(O\q a (z)b (O)\B( PB )) 



i 



j B + ^5)75 



(ki) 



{{£>b + M B ) 75 [M^i) + V2><Mki)] } 



(3a ' 



HI) 



where ki is the momentum of the light quark in the B meson, and n = (1,0, Ox), and v 
(0, 1,0t). The normalization conditions for these two wave functions are 



The relations between (f)B, <j>B and 



Sb 



J d 4 fci0g(ki) 



fs 



:il2) 



Bi YB 



arc 



>B 



B ' 



B 



'B 



(113) 



In practice it is convenient to work in the impact parameter b space rather than the transverse 
momentum space (fcj_-space). So we make a Fourier transformation J d 2 k±e~ l]iA -' h to transform 
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the wave functions and hard amplitude into 6-space. Then the normalization condition for 4>b 
and 4>b can be expressed as 



dxcpB^x, 6 = 0) 



h 



2 v / 2lVe 



dx(pB(x, b = 0) = . 



(114) 



4>B and 4>b include bound state effects; they are controlled by nonperturbative dynamics. They 
can be treated by models or by solving the equation of motion in heavy quark limit. 

In some particular models, (pB and 4>b can be selected such that the contribution of (pB is 
the next-to- leading-power A/ms |25j . In this case the contribution of 4>b can be neglected at 
leading-power. Hence, only 0b is considered in this case. We adopt the model for 0_b in the 
impact parameter b space, which is widely used in the study of B decays in the perturbative 
QCD approach [Tlj 



bsix, b) = Nbx 2 (1 — x) 2 exp 



1 fxM B Y 

2 \ uj b ) 



(115) 



where the shape parameter lob has been determined as ujb = 0.40GeV, and Nb is the normal- 
ization constant. 

In B —* />7 and ujj decays, p and oj meson can only be transversely polarized. We only need 
to consider the wave function of transversely polarized p or u meson. They are defined by 



(p/w(P,e£)|4(*)«/9(0)|0) 



1 



dxe 



ixP-z 



My^Wvix) 



My 



-ie, 



P ■ n + 

For the transverse p meson, the distribution amplitudes are given as |15j : 



f p {x) 



where t = 1 — 2x and 

lT 



f7 



V 



2^/2N~ c 



3/ P 
4^6 



x{l-x) [l + 0.2C 2 3/2 (t) 



-(l + t 2 ) + 0.24(3t 2 -l) 
+0.12(3 -30t 2 + 35t 4 
t [l + 0.93(10x 2 - lOx + 1) 



(116) 

(117) 

(118) 
(119) 



fv (v) 

9±> 



2^2N ( 



cj>y(x) 



fv d ( a ) 

-A- 



8^W c dx^ 



We use eoi23 = 1 and set the normalization condition about fa = {<fi±,g±,g± } 



as 



/ 

Jo 



dxfa (x) 



1. 



The Gegenbauer polynomial is defined by 



C 3 2 /2 (t) = ^(5t 2 -l) 



(120) 



(121) 
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B Some Functions 



The expressions for some functions are presented in this appendix. In our numerical calculation, 
we use the leading order a s formula as 



a s (n) 



2vr 



/3 ln(/i/A n 



0o 



33 - 2n f 



(122) 



and we fix the number of the flavor as nj = 4. The explicit expression for Sudakov factor s(t, b) 
is given by j^H] as follows: 



s{t,b) 
A 
B 



1 dfj, 

1/6 V 

O'.s 



hi 



A(a s {(i)) + B(a s (ii)) 



7T V 7T 



67 7T 



10 



rif + -B ln — 



27 



2« s , 
In 

3 7T 



;2 7B _ ! 



(123) 
(124) 
(125) 



where je = 0.5722 is Euler constant and Cf = 4/3 is color factor. The meson wave function 
including summation factor has energy dependence, 

(j) B (xi,bi,t) = 4> B (xi,b 1 )exp[-S B (t)}, (126) 
(j) V (x 2 ,t) = <f>v(x2)exp[-Sv(t)] t (127) 

and the total functions including Sudakov factor and ultraviolet divergences are 

S B (t) = s(x 1 P 1 ~,b 1 ) + 2 [ ^ 7 (a s (/2)), (128) 

S v (t) = S (x 2 P 2 -,6 2 ) + s((l-x 2 )P 2 -,& 2 ) + 2 I ^ 7 ( a ,(/2)). (129) 

Ji/b 2 V 

Threshold factor is expressed as below |231 127j . and we take the value c = 0.4: 

S t (x) = 2l+2 ^f /2 + c) [x(l-xW. (130) 

y/TTL (1 + C) 
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